Polarization fluctuations in vertical cavity surface emitting lasers: 
a key to the mechanism behind polarization stability 



On 
OS 

> 
o 

in 

(N 



C0 ' 

O ' 

• i-H ' 

-f— > ■ 

Or 

q 

CO 

O 

CO 
>Y 



in 
> 
o 

(N 
O 
SO 
O 
l> 
On : 

CO ' 

o : 

CO ' 
Ph. 



X 



Holger F. Hofmann and Ortwin Hess 
Institut fur Technische Physik, DLR 
Pfaffenwaldring 38-4-0, D-70569 Stuttgart, Germany 
(February 2, 2008) 



We investigate the effects of the electron-hole spin dynam- 
ics on the polarization fluctuations in the light emitted from a 
vertical cavity surface emitting laser (VCSEL). The Langevin 
equations are derived based on a rate equation model includ- 
ing birefringence, dichroism and two carrier density pools sep- 
arately coupled to right and left circular polarization. The 
results show that the carrier dynamics phase lock the polar- 
ization fluctuations to the laser mode. This is clearly seen in 
the difference between the fluctuations in ellipticity and the 
fluctuations in polarization direction. Separate measurements 
of the polarization fluctuations in ellipticity and in polariza- 
tion direction can therefore provide quantitative information 
on the non-linear contribution of the carrier dynamics to po- 
larization stability in VCSELs. 



I. INTRODUCTION 

Although vertical cavity surface emitting lasers (VC- 
SELs) are often highly symmetric around the axis of laser 
emission, practical devices usually emit linearily polar- 
ized light. The investigation of the weak anisotropics 
responsible for this polarization stability has been the ob- 
ject of a number of recent experimental and theoretical 
studies p]-0]. One of the questions raised in this search 
for a better understanding of the polarization properties 
of VCSELs is whether it is sufficient to consider the ef- 
fects of the optical anisotropics separately from the car- 
rier dynamics or whether the highly anisotropic satura- 
tion of the gain in quantum well VCSELs must be taken 
into account as well. This question is particularily com- 
plicated because semiconductor lasers are type B lasers 
and usually it is not realistic to adiabatically eliminate 
the carrier dynamics from the laser equations. The cor- 
rect description of gain saturation effects therefore re- 
quires the inclusion of carrier densities as dynamical vari- 
ables. A rate equation model for the polarization dynam- 
ics of quantum well VCSELs including the spin dynamics 
of the carriers was introduced in 1995 by San Miguel and 
coworkers [Q. This model shows that the effects of the 
linear anisotropies can be greatly modified as a result of 
the carrier dynamics A. In the case of relaxation oscil- 
lation frequencies greater than the spin relaxation rate 
the carrier dynamics can be observed as relaxation os- 
cillations in the fluctuations of ellipticity and of polar- 
ization direction Recent experimental results, how- 
ever, suggest that the spin relaxation rate is likely to be 
greater than the relaxation oscillation frequency ||(|,|]]. 



In one study, the contributions of the carrier dynamics 
have been neglected altogether since the observed polar- 
ization stability could be interpreted in terms of a type A 
laser model including only linear optical anisotropies 
A direct experimental determination of the carrier dy- 
namics contribution to polarization stability would thus 
be useful to test the validity of the different models. 

In this paper we therefore investigate the different con- 
tributions of the linear anisotropies to the polarization 
fluctuations in the case of fast spin relaxation. In section 
H the rate equations are introduced and the Langevin 
equations at the stable point are formulated. In section 
P| the Langevin equation is solved and the resulting po- 
larization fluctuations are presented. In section |^ the 
difference between the contributions of birefringence and 
dichroism and the contributions of the relaxation oscil- 
lation dynamics of the carriers is discussed and exper- 
imental possibilities of identifying the contributions are 
proposed. In section [y] the conclusions are presented. 



II. POLARIZATION DYNAMICS OF THE SPLIT 
DENSITY MODEL 

A. The rate equations 

In the model introduced by San Miguel Q] , the carrier 
density is subdivided into two carrier density pools in- 
teracting only with right or left circular polarized light, 
respectively. The physical justification for this assump- 
tion is the conservation of angular momentum around 
the axis of symmetry. In the following equations, we will 
use the parameters D for the total carrier density above 
transparency and d for the difference between the carrier 
density interacting with right circular polarized light and 
the carrier density interacting with left circular polarized 
light, n is the total number of photons in the cavity. The 
polarization is described using the normalized Stokes pa- 
rameters. In terms of the complex amplitudes of the 
circular polarized light field modes E + and EL, these 
are 



-Pi 



P?. 



P 3 = 



E*E- + E* E, 



E1E 4 



.E* + E_ 



E*_E_ 

. - E* E< 



E* + E+ + E*_E. 
E* + E+ - E*_E_ 



E*+E+ 



E*_E_ 



(la) 
(lb) 
(1c) 



1 



The relevant timescales of the laser process are given by 
the rate of spontaneous emission into the laser mode, 2w 
(usually around 10 6 — 10 7 s _1 ), the rate of emission into 
non-laser modes 7 (usually around 10 9 — 10 10 s _1 ), and 
the rate of emission from the cavity, 2k (usually around 
10 12 — 10 13 s _1 ). In addition, the spin flip scattering rate 
7 S is an important timcscalc for the polarization sensitive 
interaction between the light field and the carrier densi- 
ties. It is expected to be around 10 10 — 10 12 s _1 . In light 
of the experimental results presented in we w ^ m 

the following consider 7 S to be much faster than the re- 
laxation oscillation and adiabatically eliminate the spin 
dynamics. Note that the opposite case, i.e. where 7 S is 
much smaller than the relaxation oscillation frequency, 
has been described in ||. 

Similar to the definition of the Stokes vector, the 
anisotropies can be defined as vectors. The orientation 
of the vector indicates the polarization for which the re- 
spective physical quantities are at a maximum. There 
are three types of anisotropies: 

1. the relative gain anisotropy given by g, such that 
the rate of spontaneous emission into the laser 
mode is 2w(l + P ■ g) 

2. the loss anisotropy 1, such that the rate of photon 
emission from the cavity is given by 2k(1 + P • 1) 

3. the frequency anisotropy ft, such that the length of 

is equal to the frequency difference between the 
modes of orthogonal polarization. 

Since we only consider small anisotropies, we will neglect 
the effects of gain and loss anisotropies on the total in- 
tensity of the laser process by assuming that 1 + g • P « 1 
and 1 + g • P ~ 1. The rate equations are then given by 
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e~3 indicates the unit vector in the direction of the 3rd 
component of the stokes vector. 

\i is the injection current above transparency and a 
is the lincwidth enhancement factor which describes a 
shift in frequency due to the electron-hole density in the 
quantum well. 



laser process, we can adiabatically eliminate the carrier 
density difference d by using 
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We can also assume that d will always be much smaller 
than the total carrier density D, such that wDn + 
wdnP 3 ss wDn. The rate equations of the total car- 
rier density D and the total photon number n are then 
independent of the polarization dynamics. 
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In order to describe the polarization dynamics at con- 
stant laser intensity the stationary solutions of the to- 
tal carrier density D — 2k/ w and the stationary photon 
number n may be applied to the polarization dynamics. 
The equation for the Stokes parameter dynamics then 
reads 



dt 



P = [(s- 
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-(ft + a P 3 e 3 ) x P 



(5) 



with s = 2k (g — 1) as the total dichroism. Note that 
for 2Kwn/"f s = this equation is essentially the Stokes 
parameter version of the linear equation used in || . The 
interpretation of the polarization stability in terms of a 
type A laser model is possible because d can be adiabat- 
ically eliminated. It should be noted, however, that the 
type A model does not correctly describe the dynamics 
of the field intensity. The use of normalized Stokes pa- 
rameters is therefore more appropriate. In particular, the 
possibility of relaxation oscillations in the laser intensity 
is not included in the equations used in p. 

The effects of the split density model are given by the 
terms proportional to 2Kwn/"/ s . 2kwti is approximately 
equal to the square of the relaxation oscillation frequency. 
To estimate the importance of the split density model 
contributions to the laser dynamics it is therefore useful 
to compare the relaxation oscillation frequency with the 
spin flip scattering rate. 



B. Adiabatic elimination of the carrier density 
difference d 

If we assume that the spin flip scattering rate 7 S is 
much faster than all the other timescales involved in the 



C. Langevin equation 

In the case of stable linear polarization with parallel 
birefringence ft = Qei and dichroism s = se'i the sta- 
tionary Stokes vector is P = ei. Fluctuations of P are 
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given by the components P 2 for fluctuations in polariza- 
tion direction and P 3 for fluctuations in ellipticity. The 
linearized Langevin equation for P 2 and P 3 derived from 
equation (||) is 



dt 



Pi 
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—s — f2 — axn 
—s — \n 



Pi 
P 3 



■f(t) (6) 



where x = 2kw/j s is the coefficient of the laser intensity 
dependent contribution of the carrier dynamics in the 
split density model. 

The contributions of dichroism s, birefringence f2 and 
carrier dynamics x n can be recognized clearly in this ma- 
trix equation. The noise term f (t) is a consequence of 
the vacuum fluctuations in the electromagnetic field en- 
tering the cavity and in the dipole density of the gain 
medium. The magnitude of the noise terms for the 
Langevin equations may be derived by considering that 
photonic shot noise must be present both in the circular 
polarized modes and in the linear polarized modes at 45 
degrees to the polarization of the laser light. 



(fvs(t)f P2 (t + T)) = -S(T) 

n 



(7a) 
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The factor of An/n is the rate Ann at which photons en- 
ter and leave the cavity, divided by the squared normal- 
ization of the Stokes parameter. This is the minimum 
noise term necessary to satisfy the quantum mechanical 
uncertainty relations. Additional noise may arise from 
reabsorption of photons into the laser medium due to in- 
complete inversion. This effect will be extremely strong 
in ultra low threshold lasers p!p| . For typical VCSELs 
however, the minimal noise terms should be the main 
contribution to f(t). 



III. SOLUTION OF THE LANGEVIN EQUATION 

A. Linear response near the stationary point 

The eigenvalues A± and the left and right eigenvectors 
a± and h± of the 2x2 matrix describing the relaxation 
dynamics of polarization fluctuations in equation (||) can 
be analytically determined to obtain the linear response 
of the laser to the polarization fluctuations. 
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The eigenvalues already indicate an intensity dependent 
contribution to both the frequency and the relaxation 
rate of polarization fluctuations. The linear response to 
the polarization fluctuations is then given by the Greens 
function 
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The fluctuations in the polarization of the laser light can 
be determined by applying this Greens function to the 
noise term f(t) in the Langevin equation. 



B. Polarization fluctuations 

The fluctuations of polarization direction P 2 and ellip- 
ticity P3, as well as their correlations are given by the 
following correlation functions: 

(P2(t)P 2 (t + T))= 



(2s + x n )n 
(1 + ^K (s+x " /2)t cos(o, t) 
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The frequency luq is given by the imaginary part of the 
eigenvalues, 
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Figure [l] shows the two-time correlations as a function of 
the delay time r for a typical choice of parameters. 

In addition to the quantitative changes in the relax- 
ation rate and the oscillation frequency, the non-linear 
contribution of the carrier dynamics introduce clear qual- 
itative modifications to the polarization fluctuations. 
This is a consequence of the split density model in which 
the ellipticity is stabilized by the carrier dynamics in ad- 
dition to the dichroism. Therefore the fluctuations in el- 
lipticity are smaller than the fluctuations in polarization 
direction. Also, the phase shift between the oscillations 
of ellipticity and of polarization direction is not 7r/2 as 
one would expect for a purely linear birefringence. These 
differences between a linear optical model and the split 
density model may provide possibilities for experimental 
investigations of the split carrier density contributions to 
polarization stability. 
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C. Polarization noise in the emission spectrum 

Because the laser field amplitude E\\ is large compared 
to the noise, the two-time correlation function of the field 
amplitude polarized orthogonally to the laser mode E± 
may be determined from the two-time correlations of the 
normalized Stokes parameters. In fact, the measurement 
of the Stokes parameters corresponds to a heterodyne 
detection of the field dynamics of E± . Consequently, the 
two-time correlations of E± may be obtained using 
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The measurement of the fluctuations in ellipticity and 
in polarization direction is therefore equivalent to a phase 
sensitive measurement of the fluctuations in the laser cav- 
ity mode of orthogonal polarization to the lasing mode. 
The difference between the in-phase fluctuations cor- 
responding to the polarization direction and the out- 
of-phase fluctuations corresponding to ellipticity corre- 
sponds to a phase locking effect between the emissions 
into the non-lasing polarization mode and the lasing 
mode. 

The two-time correlation function of E± also shows 
some features of the phase locking. However, the infor- 
mation is somewhat hidden by the summation over the 
separate contributions. Neglecting the fluctuations in Eu , 
the two-time correlation function of E±_ is 
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The Intensity spectrum I±(uj) = (E*j_((j)E±((j)) is the 
Fourier transform of the two-time correlation function. 
It is given by two Lorentzians, 



2K,(£l+axn/2) 
(2s+ X n)(n+a X n) 

^ 20 
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An example of such a spectrum is given in figure us- 
ing the same parameters as for the two-time correlations 
shown in figure [l[ Note that the intensity is given in 
units of photon number inside the cavity. The intensity 
emitted is given by 2k times this value. The minimum 
noise asumption of equations ( |7a| ) and ( ffb| ) may be tested 
by comparing the total intensity in the orthogonally po- 
larized mode with the predictions of equation ( |l4|) . 

For x / the theory predicts not only the peak at 
—luq from the laser line but also a much smaller peak at 
+loq. To clarify the quantitative relation between the two 
peaks it is useful to treat Xn/f2 as a small perturbation. 
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If x n /^ is large enough the small noise peak at the op- 
posite side of the laser line should be sufficient for a deter- 
mination of x- The fact that no such peak was observed 
in indicates that x n is indeed small compared to 

il in the devices studied. 



IV. EXPERIMENTAL POSSIBILITIES 

A. Spectrum of light polarized orthogonally to the 
lasing mode 

If the carrier density dynamics is negligible for po- 
larization stability the birefringence Q is given by the 
frequency difference ujq betweeen the laser line and the 
emission line of the orthogonally polarized mode and the 
dichroism s is given by one half of the linewidth at half 
maximum Awfwhm- This is the assumption used in ||. 
That the lincwidths reported in that paper seem to be 
larger than 2s may be a consequence of the carrier den- 
sity dynamics. In particular, the linewidth with carrier 
density dynamics is given by 



Alo fwhm = 2s + xn. 



(16) 



The frequency shift ujq is much harder to identify. Ac- 
cording to the theory the frequency should increase with 
intensity until x n = 2aVL and then decrease again un- 
til overdamping occurs at x n — 2f2(a + \J o? + 1). The 
linear increase of frequency is 



luq w + axn/2. 



(17) 



This frequency shift has first been derived by van der Lem 
and Lenstra |1J] who argue that the shift is an alpha- 
enhanced saturation effect. 

The most important feature predicted by the split den- 
sity model is the peak on the opposite side of the laser 
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line. The ratio of the peak intensity at lu = ljq and at 
lu = —luq is approximately given by 

#^«®V + D. (18) 
I±(—u)q) 4S2 

If the frequency shift uiq, the linewidth IS.lufwhm and the 
intensity ratio I_l(u>o)/I_l(— u>o) are measured it is possi- 
ble to determine the contribution of the carrier density 
dynamics x n as we U as the dichroism s and the birefrin- 
gence fi at a fixed laser intensity. Note that to obtain 
correct quantitative results the linewidth enhancement 
factor a must be known. The peak at to = +luq may 
be very small because of the second order dependence on 
Xn/fl. For x n /Q = 0-1 the intensity in that peak would 
be 100 to 1000 times lower than that in the main peak 
at lu — —loq. This may be the reason why it has not 
been observed in previous experiments. If there really is 
no peak at lu = +luq however, it must be concluded that 
the split density model is not a valid description of the 
polarization dynamics in the device under consideration. 



B. Separate measurements of the fluctuations in 
ellipticity and in polarization direction 

While it is possible to observe the effects of phase lock- 
ing between the E± mode and the lasing mode En in the 
spectrum by measuring the small peak at the opposite 
side of the laser peak separate measurements of the fluc- 
tuations in ellipticity P 3 and in polarization direction P 2 
are more sensitive to the effects of \n and reveal more 
details of thephase locking. This is most clearly seen in 
figures [l] and 0. At x n /£l = 0.5 the intensity at lu = +ljq 
is still very small while the features of the polarization 
fluctuations clearly reveal strong effects of the carrier dy- 
namics. 

The most important indicator is the ratio of fluctu- 
ations in polarization direction and in ellipticity. This 
ratio can even be measured at low time resolutions since 
it is a constant over all times and frequencies. 



(P 2 (*)iM* + T)) 
(P 3 (t)P 3 (t + T)) 



= 1 



(19) 



Note that the a factor is responsible for the different 
magnitude of the fluctuations. Without the effects of 
the a factor the phase locking effects would only appear 
in the correlations between ellipticity and polarization 
direction (P 2 (t)P 3 (t + r)) and (P 3 (t)P 2 (t + r)>. If these 
correlations are measured it is interesting to determine 
the phase shift 8(f) with 
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(P 2 (t)P 3 (t + T)) + (P 3 (t)P 2 (t + T)) < 
(P2(t)P 3 (t + T))-(P 3 (t)P 2 (t + T))< 

Xn v 
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This phase shift shows how the additional damping 
of the ellipticity by the carrier density dynamics makes 
the ellipticity fluctuate less than ir/2 out of phase with 
the polarization direction. The frequency luq, the relax- 
ation rate s + x n /2, the fluctuation ratio (P 2 (t)P 2 (t + 
t)) / (P3(t)P3(t + T)) and the correlation phase 6<fi provide 
all the information needed to determine the dichroism s, 
the birefringence f2, the contribution of the carrier den- 
sity dynamics x n and the linewidth enhancement factor 
a. Note that the time resolution necessary to measure 
the dynamics of polarization fluctuations in this regime 
is given by the birefringence fi. This suggests that a 
time resolution of several nanoseconds may actually be 
sufficient. 



V. CONCLUSIONS 

The calculations presented here clearly demonstrate 
how the carrier dynamics modify the polarization fluc- 
tuations both in the spectrum of E± and in the Stokes 
parameters of ellipticity P 3 and of polarization direction 
P 2 . Even if the contribution of the carrier density dynam- 
ics x n is very small compared to the birefringence f2 a 
careful analysis of the experimental data on polarization 
fluctuations should reveal these modifications. Such an 
analysis would show whether an interpretation of polar- 
ization stability in terms of birefringence and dichroism 
only is valid or not. If the carrier dynamics of the split 
density model are relevant to the polarization stability of 
a given device, the measurement of polarization fluctua- 
tions provides information on the optical anisotropies f2 
and s, on the spin relaxation rate 7 S and on the linewidth 
enhancement factor a. If the non-linear contribution x 
turns out to be zero, it seems likely that valence bands 
with an angular momentum other than 3/2 around the 
axis perpendicular to the quantum well also contribute 
to the laser process in vertical cavity surface emitting 
lasers. 
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FIG. 1. Two-time correlations of the polarization fluctu- 
ations (Pi(t)Pj(t + t)) for a = 2 , s = 0.5 GHz, Q = 4 
GHz and \ n — 2 GHz. (a) shows the fluctuations in ellip- 
ticity (P-i(t)P-i(t + t)) and the fluctuations in polarization 
direction {P2(t)P2(t + r)). (b) shows the correlations be- 
tween ellipticity and polarization direction, (P2(t)Ps{t + r)) 
and (P s (t)P 2 (t + r)}. 



FIG. 2. Spectrum I±(u>) of the emission polarized orthogo- 
nally to the lasing mode for a = 2 , s — 0.5 GHz, Q = 4 GHz 
and x n — 2GHz. Even though x n /^ = 0.5 the peak near +5 
GHz is still very small. 
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